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Abstract 

We introduce a concept of a fractional-derivatives series and prove that any linear par- 
tial differential equation in two independent variables has a fractional-derivatives series so- 
lution with coefficients from a differentially closed field of zero characteristic. The obtained 
results are extended from a single equation to I?-modules having infinite-dimensional space 
of solutions (i. e. non-holonomic D-modules). As applications we design algorithms for 
treating first-order factors of a linear partial differential operator, in particular for finding 
all (right or left) first-order factors. 

keywords: Newton-Puiseux series for D-modules, fractional derivatives, factoring linear 
partial differential operators 

AMS classification: 35C10, 35D05, 68W30 

Introduction 

It is well-known that any polynomial equation t{x,y) = has degy{t) (counting with multi- 
plicities) zeroes being Newton-Puiseux series (see e. g. [26j ) 

yix)= Vi^'"'" (1) 

io<*<oo 

for suitable integers g > l,zo and the coefficients yi from an algebraically closed field. 

In this paper an analogue of Newton-Puiseux series for partial linear differential equations 
T = is proposed, and we prove that T = has a solution of this form. Whereas a 
Newton-Puiseux series is developed for a (plane) curve, we restrict ourselves with linear 
partial differential operators T in two derivatives dx , dy (in case of 3 or more derivatives there 
are no solutions of this form in general, see Remark I4.10p . 

One of the principal features of Newton-Puiseux series is the appearance of fractional 
exponents. Thus, a question arises, what could be an analogue of fractional powers, so to 
say "fractional derivatives"? An evident observation shows that in the derivative y'{x) = 
Ylii~'^/Q + the i-th. coefficient depends on the {i — q)-th. coefficient of y{x) itself. 
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That is why as a differential analogue of Newton-Puiseux series we suggest a fractional- 
derivatives series of the form 

0<i<oo 

where hi being elements of a differentially closed (or universal in terms of [13J) field F 
and G'(~*/'^) is called {—i/q)-th fractional derivative of G. The symbol G = G^^^ = 
G(s2,...,Si:){fiT f2, ■ ■ ■ , fk) is defined by rational numbers 1 > S2 > • • • > -Sfc > and 
/i)---)/fc & F (if to continue the analogy with curves, G plays a role of a uniformizing 
element). For any rational s the s-th fractional derivative G^*^ fulfills the identity 

dG^ = ((i/i)G(i+^) + ((i/2)G(^2+«) + . . . + (dfk)G'^'^+''> 

where either a derivative d = dx ox d = dy. The common denominator g of S2, • • • , plays a 
role similar to one of the common denominator of the exponents in a Newton-Puiseux series 
([T]). The inequality k < q holds. 

In a particular case k = 1 we have q = 1 and as G one can take g{fi) for any univari- 
ate ("undetermined") function g, provided that the composition makes sense, the fractional 
derivatives G^**-* = g^^\fi) for integers s. We note that finite sums 

E /.^G(-) 

io<«<ii 

(so, ioi k = q = 1) appear in the Laplace method as solutions of some second-order equations 
T = (see e. g. [lES]). 

One can find necessary in the sequel information on D-modules in [2], [IS], a survey on 
their algorithmical aspects in [20]. We mention that there are also applications of Newton 
polygons over the Weyl algebra C[x, dx]'- in [16| to meromorphic connections, in [IT] to micro- 
differential operators and in [18] to the Fourier transform. In case of linear ordinary differential 
operators Newton polygons are employed to produce the canonical form basis of the space of 
solutions (see e. g. p7], also ffj where an algorithm for this problem with a better complexity 
bound was designed). A similar form of solutions for linear partial differential operators were 
studied in [l] where, nevertheless, also examples are exhibited of operators without solutions 
of this form. On the problem of factoring a linear ordinary differential operator one can look 
in [19], see also [7j. 

In Section [1] we introduce the principal concept of fractional-derivatives series and give 
some their basic properties. 

In the sequel the crucial role plays the multiplicity m of a linear factor of the symbol of 
the linear partial differential operator T (with coefficients in F) of an order n: the symbol is 
a homogeneous polynomial in two variables dxfi,dyfi of the degree n which corresponds to 
the highest derivatives of T. In Section [2] we develop a method for constructing fractional- 
derivatives solutions of r = and prove the existence of such a solution with q < m. The 
method is similar to the Newton-Puiseux expansion, it produces a relevant convex polygon 
similar to the Newton one, but differs in several aspects. The main of the latter is that the 
leading equation corresponding to a certain (leading) edge of the polygon is not a univariate 
polynomial unlike the Newton-Puiseux expansion, but rather a non-linear first-order partial 
differential equation. This creates difficulties in defining a multiplicity of a solution of the 
leading equation. Also it is unclear, what could be a differential analogue of the statement (cf. 
above) that an algebraic equation t = has precisely degy{t) Newton-Puiseux series solutions 
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of the form ([T])? Partially these questions are answered for the introduced in Section [3] generic 
fractional-derivatives series solutions. 

In Section S] the result of Section [2] is extended from a single partial linear differential 
equation to a system of equations in several unknown functions having an infinite-dimensional 
space of solutions (or in other words, to a D-module of a non-zero differential type, one can call 
it a non-holonomic I?-module). To this end for any left ideal J C F[dx,dy\ of the differential 
type 1 we yield an operator p £ F[dx,dy\ and show that any fractional-derivatives series 
solution of the equation p = which corresponds to a linear factor (different from dyfi) of 
the symbol of p, is a solution of the ideal J as well. In Section [5] we exploit the relation of 
equivalence of ideals introduced in [ID] and establish a kind of duality between equivalence 
classes of non-holonomic ideals and their sets of fractional-derivatives series solutions. Namely, 
it is proved that two non-holonomic left ideals J, Ji C F[dx,dy\ are equivalent if and only 
if their respective sets of fractional-derivatives series solutions coincide. Also we express the 
quotient of the spaces of fractional-derivatives series solutions of non-holonomic ideals J C Ji 
via the module of relative syzygies [TU] of this pair of ideals. 

In Section [6] it is shown that in case of a separable operator T any its power series 
solution can be obtained as a sum of specifications of its suitable fractional-derivatives series 
solutions, thereby establishing completeness of the latter. In Section[7]we provide applications 
of fractional-derivatives series to studying first-order factors of an operator, exploiting that 
in case of a first-order operator T = dy + adx + h its fractional-derivatives series solutions 
turn to a single term of the form hG{f) where {dy + adx)f = and h being a particular 
solution of r = 0. In Subsection 17.11 an algorithm is designed which finds first-order factors 
of a given operator, and in Subsection 17.31 an algorithm which constructs the intersection 
of all the principal ideals generated by the first-order factors of the operator. In Section [H] 
the possible fractional-derivatives series solutions of a second-order operator obtained by the 
algorithm from Section [5] are described. This description can help to imagine the shape of 
fractional-derivatives series solutions and the difficulties which appear while their developing. 

1 Fractional-derivatives series 

Let be a differential field of the characteristic with the derivatives {dj} and a subfield of 
constants C C F |13j . 

Definition 1.1 Let /i, . . . , /^^ E F and rational numbers 1 > S2 > • • • Sfco > 0. We introduce 
a symbol G = G^^^ = Gs2,...,SfeQ (/i, f2, ■ ■ ■ , fko) together a set {G^^^}s^q of its fractional s-th 
derivatives satisfying the following rule of differentiation for any derivative d = dj : 

dG'-'^ = (d/i)G(i+") + (4f2)G("2+") + • • • + (d/fco)G(^^o+^) 

Clearly, these differentiations commute with each other and one can consider the free 
F-module with the basis {G^'^^jsgQ as a Z)-module. 

Definition 1.2 Let q be the common denominator of S2, ■ ■ ■ , s^g and hi £ F,i > 0, sqq £ Z. 
Then 

H= hiG^'"-'/''^ (2) 

0<i<oo 

we call a fractional-derivatives series. 
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For a given G all the fractional-derivatives series (with added 0) constitute a L'-module 
(we study it below in Section H]). Obviously, ko < q. 

It is easy to see that G satisfies a suitable linear partial differential equation with coeffi- 
cients in F. 

Remark 1.3 The symbol G plays a role in H similar to the role of the parameter x in a 
Newton- Puis eux series (Cp. In particular, specifying the values of x in a certain field one gets 
points of (a branch of) the curve given by ([2]j. Here one can also provide some specifications 
of G. Indeed, for an arbitrary family {ci/q}i£z where Ci/q € C the following set 

_ \^ „ £1 lko_ 

^ - c-s-n-ns2 jfcQSfeQ . , . , 



satisfies Definition \l.l\ 

For example, in case when F is the ring of functions analytic in a certain neighborhood of 
a given point in the multidimensional complex space and the absolute values |cj/^| are bounded, 
the latter series also converges in a suitable neighborhood. 

From now on let F have two derivatives {dx, dy}. Consider a linear operator 

T = To + • ■ ■ + r„ (3) 

of the order n where Tp = X]o<j<p bj^pdidy~-^ contains the derivatives of the order p and 
the coefficients bj^p G F. The following lemma holds, in fact, for an arbitrary number of 
derivatives, nevertheless, the assumption that F has two derivatives simplifies the notations 
and in the sequel we deal just with operators in two derivatives (one can verify lemma by a 
direct calculation). 

Lemma 1.4 dxdy~-' (hG) equals the sum of the terms of the form 



J \ P-J 



Wi \ ■ • • Wj\[\ , • • • , ll,mi 5 • • • ) ^fco,l) • • • 5 ^A;o,m(;y ) v V^l,li • • • ) 5 • • • i ^fco,l ) • • • ) ^fccm^jj j ^ J 

l<i<mi l^^^Tifcg 

for all partitions H h li,mi H h lko,i H \- ho,mk^ + I = j of j and ri,i H h 

ri,mi + • • • + rko,i + • • • + rko,mkg + r = p- jofp-j such that l^^i + r^^i > 1 for ev- 
ery 1 < K < ko, 1 < i < tUk , where wi, . . . denote the cardinalities of the partition 
of the triples (/i,i, n^i, 1),. . . ,{li 

,1) ''"fco,!) ^0)5 • • • ) {^0,171^ ) 1^ko,mk,. ) ^0) 

into equal ones, in particular, wi + ■ ■ ■ + wn = nii + • • • + m^Q. 



2 Constructing fractional-derivatives series solutions 

Prom now on we suppose that the field F is differentially closed (or universal in terms of ) . 

The main purpose of this section is to prove that a linear partial differential equation 
T = 0, see has a solution of the form To simplify the notations we put sq = and 
h = ho^Oml^. 
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Denote by Tp{dxfi,dyf2) = YjQ<j<p bj,p{dxfiy{dyfiy ^ a homogeneous form of the degree 
p in dxfi,dyfi. Sometimes, r„ = symb{T) is called the symbol of T. Fix a linear factor 
o-idxh + CL2dyfi of Tn having a multiplicity m, the coefficients ai,a2 G F. 

Expanding T{H) with respect to the fractional derivatives {G^^^}s for A: = 1 (in other 
words, assuming for the time being that dG^'^ = {dfi)G^^+''\ see Definition [H]), we get 
that the coefficient at G^") vanishes, i. e. h ■ symb{T) = 0. Thus, we can suppose that 
{aidx + a2dy)fi = 0. Choose any such /i with grad{fi) / 0. 

For k > 2 we introduce an auxiliary polygon playing the role similar to the Newton 
polygon. Now let /c = 2, in other words, we assume (for the time being) that dG^"^ = 
(d/i)G(i+") + {df2)G^'^+'\ The next purpose is to construct S2 and f2- It suffices to consider 
the expansion of the first term T{hG) of T{H) (we'll come back to this issue at the end of the 
present section). When we talk about the expansion of T{hG) we always refer to Lemma [1.41 
If a term &(ni<i<t('^a-'^j^'/2))G'^*^*^*^ occurs in T{hG), where 6 is a differential polynomial in 
/i and in h (being linear in h), then we place the point (s, t) in P2. As P2 we take the convex 
hull of these points with the origin (0, 0). If to assign the weight 1 to every derivative d^^dyfi 
then any term in b gets the weight s due to Lemma 11.41 

One can observe that P2 lies to the left from the line Li = {s + t = n} with the slope 1 
(under the slope of the line {s + jt = const} we mean j) again due to Lemma 11.41 Moreover, 
the point (n — m, m) € Li belongs to P2 because the non-zero term 



, w,^ • + a2dy)f2r ■ 

((ai4 + a2rfs/)/i)™ 



■m+S2m) 



occurs in the expansion of T{hG), taking into account that the factor {aid^ + a2dy)fi has 
the multiplicity m in T„, and no other term from this expansion gives a contribution in the 
coefficient at the point (n — m,m). Similarly, one verifies that the points (n — t,t) with 
0<t<m — Ido not belong to P2 ■ 

Now we assign a (yet unknown) weight S2 to every derivative d^^dyf2- Therefore, to find 
S2 < 1 we consider the edges of P2 with the positive slopes less than 1. Choose any such edge 
L2 (we call it leading) with the endpoints (ji,ii), (^2)^2); ti > t2] we have seen already that 
ti < m. Then the slope of L2 provides S2 = (ji — 32) /{ti — ^2)- 

To find /2 we consider the leading differential polynomial (52 (/2) which equals the sum 
of the coefficients at all the points of P2 which lie on L2. Then Q2(/2) coincides with the 
coefficient at G^^'i+'^^ti) the expansion of T{hG). As /2 € -F we take a solution of the leading 
equation (52 (/2) = 0. Evidently, ji + S2ti < n since the point of intersection of the line L2 
(which contains the edge L2) with j-axis {t = 0} is located to the left of the intersection of 
Li with J -axis. 

Thus, we are able to formulate the recursive hypothesis of the procedure under description 
which constructs 1 > S2 > S3 > • ■ ■ and /i, /2, /s, . . .. Suppose that S2, ■ ■ ■ ,Sk and /i, /2, • • • , /fc 
are already constructed. In addition, a polygon P^, is constructed being a convex hull of the 
points (j, t) (together with the origin (0,0)) such that a term 



"(n 



l<i<t 



4»<»/fc)G(^'+^^-*) (4) 



occurs in the expansion of T{hG) under the assumption dG = {dfi)G^^^ + {df2)G^^'^^ + • • • + 
{dfk)G^^'^\ see Definition ll.l[ A certain leading edge L^. of P^ is chosen with a slope Sk > 
and with the endpoints (^'3,^3), {ji^t/i), > t^. We name (j3,t3) the pivot of Lk and t^ the 
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multiplicity of L^. The leading differential polynomial Qkifk) equals the sum of the coefficients 
at all the points of which lie on L^.. Then Qkifk) coincides with the coefficient at G^^^'^^''^^^ 
in the expansion of T{hG). As G F a solution of the leading equation Qkifk) = is taken. 
The points of intersections of the lines Li,L2,... with j-axis decrease. Denote by the 
common denominator of S2, • • • , Sfc, obviously qi = 1. 

To carry out the recursive step, we make the assumption dG = idfi)G^^^ + idf2)G^'^^^ + 
■ ■ ■ + idfk)G^^''^ + idfk+i)G^'^''+'^\ The boundary of the polygon P^+i above the pivot of Lk 
(including the pivot itself) is the same as of P^- 

Let us calculate the points of Pk+i located on the line Lk- Denote by Bt it^ < t < t^) the 
coefficient of at the point (i3 + Sfc(t3 — t), t) G Lf^. Then Qfc = Ylt3<t<t4 ^t- observe 
that Bt contains no higher derivative d^r^dyfk with l+r > 2. Indeed, if otherwise Bt contained a 
term of the form then the coefficient of Pk at the point (j3+Sfe(i3— i),i+X^i<j<j(^i+?'j — 1)) 
would contain the term 



due to Lemma [L^ hence the point (j3 + Skits — t),t + X^i<j<t(^j + ri — 1)) should belong to 
Pk which leads to a contradiction when Yli<i<ti^i + ri — 1) > 1. 

Besides, Bt is a linear form in the derivatives of h. We claim that Bt = hBt for an appropri- 
ate differential polynomial Pt in /i, . . . , /fc. Indeed, if Bt contained a term id^^dyh)hG'^^^'^'^^^^^ 
with l+r >1 for a certain b being a differential polynomial in /i, . . . , Z^, then the coefficient of 
Pk at the point ij?, + Skit^-t),t+l+r) would contain the term /i5(4/fc)'(dy/fc)''G(^'2+'*''(*3+'+'')) 
due to Lemma [L31 therefore, the point (js + Sfc(*3 — t),t + I + r) should belong to Pk, the 
achieved contradiction proves the claim. 

Thus, Bt will be treated as a homogeneous (of the degree t) polynomial in d^fk-, dyfk- For 
more generality of the auxiliary results below we deem that Bt is a homogeneous polynomial in 
the variables vi, . . . , Vp, thereby p = 2 and vi = dxfk, V2 = dyfk- We denote the corresponding 
derivatives vi = m = dyfk+i- 

Remark 2.1 Since the main purpose of the present section is to prove the existence of so- 
lutions of the form (0j of an equation T = (see ^) it suffices to study only the canonical 
solutions, namely, when each Sk is the slope of a certain edge of Pk and fk satisfies a leading 
equation. Alternatively, one could take Sk to be the slope of some line passing through a single 
vertex, say (^'3,^3) of Pk- In this case Bt.^ifk) = 0, because Bt.^ is a homogeneous polynomial 
in dxfk, dyfk, we get that fk fulfills a certain first-order linear equation bidxfk + b2dyfk = 0. 
There is no way to bound the denominators Sk for non-canonical solutions the number of 
steps ko, moreover, the procedure of constructing 1 > S2 > S3 > • • • and fi, /2, fs, - - - could last 
infinitely- One might even choose real exponents Sk (cf. 111^ where an analogue of Newt on- 
Puiseux series solutions with real exponents was studied for non-linear ordinary differential 
equations). 

Denote by Bti^ 1^ t < t^) the coefficient at the point (js + Skit^ — t),t) G Lk of Pk+i- 
Taking into account the assumption on dG and Lemma 11.41 we have 



Kdxfk) 





Ip — I 
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Therefore, Bt = hBt where Bt can be treated as a homogeneous polynomial in vi, . . . ,Vp of 
the degree t with the coefficients being differential polynomials in fi, ■ ■ ■ , fk- Let to be the 
minimal t such that Bt / 0. Then to > 1 because Qkifk) = 0, and to < *3 because Bt^ is 
obtained from Bp^ by means of replacing Vi for Vi,l < i < p. One can view to as a kind of 
multiplicity of the solution fk in Qk- 



Lemma 2.2 to < t4 + 



(t3-f4)gfc- 



1 



<?fc 



Proof. Suppose the contrary. First we observe that the gap between the ordinates of any 
pair of consecutive points on is at least q^/qk-i and that e = (ts — t/C}qk-i/<lk is an integer 
(cf. |26j). Hence Lk contains at most e + 1 points. Without loss of generality for the sake of 
conveniency of notations we assume that Lk contains exactly e + 1 points (some among them, 
perhaps, with zero coefficients Bt). 

Due to the supposition and the choice of to we have Bt = for t4 < t < t4 + e, i. e. all 
the derivatives 

d'Qk 
dv\^ ■ ■ ■ dvi" 

of the order t vanish. Fix for the time being non-negative integers ji, ■ ■ ■ ,jp with the sum 
Ji H V 3p = ti. Then 



ih - jiY ■ ■ ■ (ip - jpV- dv\' ■ ■ ■ dvp" 



0= > " . v\'-''---Vp 

h>ji,--;ip>jp;h+---+ip=t 



^{l-U)\ d'^Bi 



l>t 



(/-t)! dvi'-'-dvi" 



due to the Euler's formula. The latter equalities can be treated as a linear (e + 1) x (e + 1) 
system with a non-singular matrix. Its non-singularity is justified by the following result [15j: 
if rii > ■ ■ ■ > rir > 0; mi > • • • > nir > 0; ni > mi, . . . , n,. > rrir then the r x r matrix with 
the entries (^',) is non-singular. Therefore, 

d^^Bi 



dvi' ■ ■ ■ dvjr 

for any / and any ji, . . . ,jp with Ji + • • • + Jp = t4, in particular Bt^ vanishes identically, the 
obtained contradiction proves the lemma. ■ 

Corollary 2.3 // to = ta then the denominator Qk-i = Qk does not change. 

Now we are in position to continue the recursive step of the procedure constructing 
Sk+i,fk+i- The polygon Pk+i either contains the edge with the slope Sk and with the ordi- 
nates to < ts, respectively, of its endpoints, or the edge of Pk+i with its above endpoint (j'a, ts) 
has the slope less than Sk. In the first case as a leading edge Lk+i one takes an edge of Pk+i 
having a positive slope Sk+i with the ordinate ts of its upper endpoint (j^jt^) less or equal 
to to- In this case (js, ts) plays the role of a new pivot with ts being the multiplicity of Lk+i. 
As above one produces the leading differential polynomial Qk+i{fk+i) and as fk+i chooses a 
solution of the equation Qk+i{fk+i) = 0. In the second case the denominator qk = qk-i does 
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not increase due to Corollary 12.31 and as Lfc+i one takes an edge of Pk+i having a positive 
slope Sfc+i with the ordinate of its upper endpoint (the pivot) (j^jt^) less or equal to t^. 
The rest is similar to the first case. 

Thus, we have described a recursive procedure constructing 1 > S2 > S3 > • • • and 
/i) 72, /sj • • • which one can view as a tree. 

Lemma 2.4 i) The common denominator q of S2, S3, . . . does not exceed 2™-"^; 

ii) there exists a branch of the tree in which the common denominator q is less or equal 
to m; 

Hi) every branch of the tree after at most of q steps arrives to a leading edge with a 
non-positive slope. 

Proof. First we recall that the multiplicity of any leading edge in P2 is less or equal to 
m. Therefore, i) follows from Lemma l2.2t if at a certain step the common denominator qk-i 
is multiplied by q^/qk-i then the multiplicity decreases at least by qk/qu-i — 1- After the 
multiplicity reaches 1, the denominator does not change anymore. 

ii) Let us take at each step of the described recursive procedure the leading edge with the 
least possible slope, while the latter is positive. The ordinate of the lower endpoint of this 
edge t4 = 0. Therefore, Lemma 12.21 entails that to ^ t^Qk-i/ Qk-, this implies ii). 

iii) follows from Definition 11.11 because < q. n 

Assume now that Pt+i in the described procedure contains an edge having a non-positive 
slope (see Lemma [23] ii) ) . Take such edge L = L^+i with the largest possible non-positive 
slope in Pk+i- We have shown above that the coefficient Bt^ at the pivot (j5, ts) of L^+i equals 
to hI3 where 13 is a suitable homogeneous polynomial of the degree ts in dxfk+i, dyfk+i with 
the coefficients being differential polynomials in /i, . . . , /fc. Denote by B the coefficient at the 
point (j5, 0) of Pk-\-i, being a linear homogeneous operator in h (one can show that the order 
of B does not exceed ^5 in the same manner as it was shown that Btr^ has the order in h). 
If B contains a term h{dxfk+iY{dyfk+iY^~^ for some / and b being a differential polynomial 
in /i, . . . , /fc, then B contains the term b{{dxy {dyY^~'- h) due to Lemma \TM Hence the order 
of B is greater or equal to t^ (actually, equals t^ as we have seen, although we use below only 
that the order of B is positive). In particular, the slope of L equals 0, and Pk+i contains no 
edges with negative slopes. In the construction under description fk+i does not appear and 
as h £ F we take a solution of the linear homogeneous differential equation B{h) =0 (which 
can be viewed clS cl leading equation on h). 

This completes the construction of the first summand hG of the solution H of the form 
([2]) . To obtain the next coefficient hi of H we observe that in the expansion ofr(/iiG(-i/9)) m 
the fractional derivatives {G^*/'^''}_oo<i<oo the highest non-zero term equals B{hi)G^^-'~^/'^\ 
taking into account that this expansion is obtained by means of the shift by — l/g of the 
expansion of T{hG) while replacing h for hi. Therefore, for hi £ F we get a linear partial 
differential equation (not necessary, homogeneous) of the form B{hi) = f (so, of the same 
order ^5) for an appropriate f € F being a differential polynomial in h, fi, . . . , fk^ (in the 
above notations ko = k). In a similar way one obtains consecutively /i2, /13, . . .. 

Summarizing, the following theorem is proved. 

Theorem 2.5 Any linear partial differential equation T = of an order n (see ^) for each 
linear factor {aidxfi + a2dyfi) of a multiplicity m of its symbol symb(T) has a non-zero 
fractional- derivatives series solution of the form with the denominator q <m. 
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One can continue every branch of the tree of the described procedure constructing 1 > S2 > 
S3 > • • • and /i, /2, • • • to a solution of the form ^ of T = 0, and every solution of the form 
^ constructed by a described procedure has the denominator q < 2"^~^ . 

Corollary 2.6 // an LPDO T has no fractional- derivatives series solutions of the form ^ 
corresponding to a factor aid^fi + a2dyfi of a multiplicity m of the symbol symb{T) with the 
denominator q < m then T is irreducible in F[dx,dy]. In particular, if T of an order n has 
no fractional- derivatives series solutions with the denominator q < n then T is irreducible. 

Remark 2.7 The bound q < m is sharp as shows the following example. Take T (see ^B^) 
such that {aidxfi + a2dyfi) has the multiplicity m in Tn, the multiplicity greater or equal to 
m — i in Tn-i for every 1 < i < m — 2 and the multiplicity in Tn-m+iji^Gspectively. Then 
the polygon P2 has the edge with the endpoints (n — m,m) and (n — m + 1,0) which being 
taken as a leading one (actually, there is no other choice for a leading edge), provides the 
slope S2 = 1/m. 

Remark 2.8 Theorem \2.5\ states the bound q < m for a particular solution. It is unclear how 
sharp is the bound q < 2"^"^ for all constructed solutions. The natural question is whether 
one can improve it by m ( one can verify it for m < 7 by the direct calculations ) ? This would 
be similar to the algebraic situation in which such a bound on the common denominator in all 
Puiseux series ([ip is well known (see e. g. 126]). We also mention that for solutions in the 
canonical form basis l2T^ of linear ordinary differential equations a similar to the algebraic 
situation bound on the common denominator (of the rational exponents) was established in 

3 Multiplicity of generic fractional- derivatives series solutions 

In the described recursive construction was chosen as a solution of the equation Qkifk) = 
0. Different choices of fk could yield different polygons Pk+i- Therefore, the set of (even 
canonical fractional-derivatives series, see Remark 12. ip solutions of the equation T = is 
quite vast. An interesting open question is whether it is possible to introduce a concept of a 
multiplicity of a set of fractional-derivatives series solutions and relate it to m? In the present 
section we give a partial answer to this question for the so-called generic solutions. 

We view Qk as a polynomial in two variables vi = dxfk,V2 = dyfk- Note that this 
polynomial is not homogeneous, consider its factorization = ■ ■ ■ P^^' P over F where /3 
is homogeneous and Pi, . . . , f3i are irreducible non-homogeneous. In the recursive construction 
from Section [2] we distinguish a case which we call generic, namely, when Pi{dxfk,dyfk) = 
for a certain 1 < i < I such that rui = minjmi, . . . , m/}, and the point {dxfk,dyfk) is a 
non-singular one of the plane curve Qk = 0. In the generic case for the multiplicity of fk we 
have to = rui due to ([5]). One can assign the multiplicity to to the set of all fk satisfying the 
generic case. We call a solution ^ generic if for each of /2, . . . , fko the generic case happens 
in the construction of ([2]). When = 1 we call ([2]) generic as well. At the end of developing 
any generic solution we arrive to a polygon Pk+i having a leading edge Lfc+i with the slope 0. 
Let the upper endpoint (pivot) of Lk+i be {j^,t^), then to this generic solution we assign the 
multiplicity t^. Observe that we have assigned the multiplicity to the set of all the generic 
solutions ^ which follow the same branch in the tree of the construction from Section [5J 



9 



Proposition 3.1 Any linear partial differential equation T = of an order n 

i) has a generic solution of the form 

ii) the sum of multiplicities of the generic solutions does not exceed n; 
Hi) the denominator of every generic solution is less than n'^(^°s'>^) ^ 

Proof. Each Q/., k >2 is non-homogeneous, that is why i) is justified taking into account 
Theorem 12.51 

ii) follows (similar to the algebraic Newton-Puiseux series [26]) by inverse induction along 
the tree of the procedure described in Section [2] due to the inequality mi + • • • + < ts — ^4. 

The latter inequality together with Lemma 12.21 imply that to ^ ^^Ig^T^g^TT' therefore, in 
developing a generic solution by means of the procedure from Section [2] there are at most 
I0S3/2 steps at which the denominator augments. At each such step the denominator grows 
less than in n times (cf. the proof of Lemma l2. 20 . this entails iii). ■ 

Remark 3.2 In a particular case m = 1 we have q = Uq = \, all the solutions of the form |^ 
are canonical, the polygon P2 contains a single edge with a slope less than 1, namely, the edge 
with the endpoints (n — 1, 1) and {n — 1,0) having the slope 0. It provides the leading linear 
equation on h of the first order, the leading equation (aidx + a2dy)fi = on fi is linear and 
of the first order as well, thus, the multiplicity 1 is assigned to the set of (generic) solutions 
in case m = 1. 

Remark 3.3 While in (0) we consider series with decreasing orders of derivatives ofG, one 
can easily verify that an equation T = for an arbitrary f (z F has a solution of the form 

0<i<oo 

where G = G{f), with increasing orders of derivatives ofG. Thus, continuing the analogy with 
plane curves, the latter series could he viewed as corresponding to expanding at finite points 
all being regular (so, without proper fractional derivatives, i. e. ko = 1 in Definition 
while (0) corresponds to expanding at the infinity. 

4 Fractional-derivatives series solutions of non-holonomic 
/^-modules 

First let J = {pi, . . . ,pi) C F[dx, dy] be a differential (non-holonomic) left ideal of the differ- 
ential type 1 |13[ I14j. This means that the Hilbert-Kolchin polynomial Kj{z) = ez -\- cq of 
J has the degree 1. Denote by symh{J) C F[dxfi,dyfi] a homogeneous ideal generated by 
the symbols of elements of J (cf. Section [2]). Then Kj coincides with the Hilbert polynomial 
Ksymh(j) El EI] (one can also deduce this from the Janet base of J [22l [10], we mention that 
the concept of Janet bases was a differential historical predecessor of the one of Groebner 
bases). Denote g = GCD{symb{J)) G F[dxfi,dyfi]. 

Lemma 4.1 The degree e of the ideal symh{J) coincides with deg{g). 

Proof. Since symb{J) C {g) it suffices to verify that dimp{{g) / symb{J)) < 00. Nullstel- 
lensatz entails that {symb{J) / g) D {dxfi,dyfi)^ for a suitable s, therefore, the homogeneous 
component 

{9)deg{g)+s = 9 ' {dxfi,dyfiY C g ■ {symb{J)/g) = symb{J) m 
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Remark 4.2 // ideal J C F[dx,dy] is holonomic (so, with differential type 0) then 
GCD{symh{J)) = 1. 

The degree e (being the leading coefficient of the Hilbert-Kolchin polynomial) is called 
the typical differential dimension of J [ T3| [Hj . 

For any homogeneous polynomial € F[dxfi,dyfi] and a € F denote by multa{go) the 
multiplicity of the linear form dxfi +adyfi in qq. Also for any p E F[dxfi, dyfi] we denote for 
brevity multa{p) = multa{symh{p)). W.l.o.g. assume that dyfi does not divide g (otherwise, 
one can perform a suitable C-linear transformation of dx,dy). We have multa{symb{J)) = 
multa{g) and e = ?7iu/ta(5) (cf- Lemma l4.ip . For the time being fix a G F such that 
multa{g) > 1. 

Now we introduce the ring R = F[dx, dy\{F[dy\)^'^ of partial-fractional differential opera- 
tors |H]. Its elements has the form pob'"^ where po G F[dx.: dy], 6 S F[dy\. One can verify (see 
[8]) that R is an Ore ring [2], any element of R can be written in a form for appropriate 
p £ F[dx,dy],b e F[dy]. Thereby R = {F[dy])~'^F[dx,dy] and pob~^ = b'^p if and only if 
bpQ = pb. Also in [8] one can find the algorithms for addition and multiplication of elements in 
R. Any element from R can be written in the form b^^ Y^(}<i<w bid\. for suitable 6, bi G F[dy\ 
(because a finite family of elements from R has a common denominator which belongs to 
F[dy\, see [8]). 

For the time being fix G = ^(^2 ,,,^s^)(/i, . . . , fk) such that {dx + ady)fi = 0, dyfi ^ (cf. 
Section [1]). Denote by F = Vg the F[dx, dy]-module which consists of all fractional-derivatives 
series of the form ([2|) added by 0. 

Lemma 4.3 V is an R-module 

Proof. For any 0^ H eV and / 6 G we claim that bH ^ 0. Indeed, let 

H = hG^'^ +Y,hiG^"''^'^\h^O;b = tnd'^+ ^ ti4,i„y^0, 

i>l 0<j<n-l 

then 

bH = /it„(dJi)"G["+") + / 

i>l 

For any Hi G y we need to prove the existence oi 11 € V such that b^^Hi = H, i. e. Hi = 

bH (the claim above implies that H is unique). Let Hi = /ii^oC^^^+Zii^iG^^^^/'^^H ; Hi^q ^ 0. 

Then we look for H = hG^^~'^'> + hiG^^~^~^/'^'> + ■ ■ ■■ Comparing the coefficients of Hi and 
bH at G^'^), we get hifl = htn{dyfi)'^ which yields h. Comparing the coefficients at G^'^ ^^'^^ 
yields hi and so on. ■ 

Remark 4.4 By the same token multiplying 6?/ 7^ p G F[dx,dy] on V is an isomorphism, 
provided that {dx + ady)fi does not divide symb{p). 

The ring R is left-euclidean (as well as right-euclidean) with respect to dx over the skew- 
field F[dy]{F[dy])~-^ , cf. Lemma 1.3 [8J. Hence the ideal J = (pi, . . . ,pi) C R is principal, let 
J = {p) for an appropriate p G J C F[dx,dy]. Then for any pQ £ J (actually, moreover for 
Po G J) the equalities 

PoP = boPo,bp= Y PjPj (6) 
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hold for suitable pj G F[dx, dy]; 7^ b, bo G F[dy]- 

According to © we have symb{po)symb{p) = symb{bo)symb{po), whence 

multa{po) = multa{bo) + multa{po) = multa{po) + multaip) > multa{p), 

therefore, multaip) — rnulta{J) since 

multa{J) = multa{g) = mmmulta{po)- 

On the other hand, from ^ we get 

multaip) = multaib) + multaip) = multai PjPj) ^ multai-J) 

Thus, the following lemma is proved. 

Lemma 4.5 For any a G F we have multaip) — multaiJ). 

Proposition 4.6 A fractional- derivatives series H with dyfi 7^ (see is a solution of 
the linear partial differential equation p = if and only if H is a solution of the ideal J. 

Proof. If pH = then from © we have = popH = bop^H. Hence p^H = due to 
Lemma 14.31 The inverse statement follows again from ([6]) and Lemma 14.31 ■ 

Corollary 4.7 For any 01,02 G F such that mult(^aida:fi+a2dyfi)isy^biJ)) > 1, the ideal J C 
F[dx,dy] has a solution of the form with a denominator q < mult(^aidxfi+a2dyfi)i^y^biJ)) 
and aidxfi + a2dyfi = 0,gradifi) / 0. 

Proof. It follows from Lemma 14.51 Proposition 14.61 Lemma 12.41 and Theorem 12.51 ■ 

Remark 4.8 If for every a (z F the ideal J has the multiplicity multaiJ) < 1 then all the 
solutions of J of the form ^ are canonical, and J has precisely e = ^a&F ''^'^^^aiJ) (which 
equals the typical differential dimension of J, cf. Lemma \4^ families of fractional- derivatives 
series solutions. Moreover, to each of these families a multiplicity 1 can he naturally assigned 
(cf. Remark\3^. 

Finally, let U C iF[dx,dy\y be a (non-holonomic) F[dx, dyj-module of the differential type 
at least 1, obviously, the differential type does not exceed 2 (recall that the differential type 
equals the degree of the Hilbert-Kolchin polynomial of U [131 114j). Denote hy ui, ... ,ui a 
free base of iF[dx,dy]f . For any 1 < r < Z consider the submodule Ur = {'}Zr<i<iPi'^i ^ ^} 
where pi G F[dx,dy\. Denote by Jr = {pr} C F[dx., dy] the left ideal being the projection of Ur 
on the r-th component. Then the differential type of U coincides with the maximum of the 
differential types of {Jr}i<r</ (one can verify this, e. g. using the Janet bases of {Jr}i<r</ 
which provide a triangular Janet base of U). Take the minimal tq such that J,.„ has the 
differential type at least 1. 

One has the natural action U y^V^ ^ V on the free F[dx, dyj-module = {X]i<i<« FliVi} 
where Hi G V (cf. Lemma HTS]) and vi, . . . ,vi is a free base of If "^KiKiPiFIi = then 
we call Yli<i<i HiVi a solution of Yli<i<iPi'^i i^^ shall choose G and thereby, V = Vq later). 
We are looking for a solution of the form Y2i<i<i FtiVi of the module U. 
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First we put Hrg+i = ■ ■ ■ = Hi = and as Hr^ ^ take a fractional-derivatives series 
being a solution of the ideal J^-q according to Corollary 14.71 in case when the differential type 
of equals 1. When the differential type of J^q equals 2, in other words, J^q = 0, we take 
as 7^ an arbitrary fractional-derivatives series. In both cases H^^^Vr^ is a solution of the 
submodule Ur^- Thus, we have chosen G = G(^s2,...,Sk)ifi^ ■ ■ ■ ^ fk) ^-nd thereby, V = Vq. As 
above we can assume w.l.o.g. that in the equation {a\dx -|- a2dy)fi = we have oi / 0, so 
{dx + ady)fi = (performing if necessary a suitable C-linear transformation of dx, dy). 

Now we construct Hj. by recursion on tq— r > 0. Suppose that we have already constructed 
an element 'Ylr+i<i<i ^i'^i being a solution of Ur+i for some r + 1 < vq. Since Jr has the differ- 
ential type (due to the choice of ro), Jr contains a certain element 7^ & € F[dy]. Consider 
a corresponding element u = bu^ + Yl,r+i<i<iPi'^i ^ ^r- According to Lemma ll?3l one can find 
Hj. such that hHr + J2r+i<i<iPi^i ~ ^- -^o^ element u = Ylir<i<iPi'^i ^ applying 
the left euclidean division in R one can represent u = prh~^u + u ioi an appropriate u G f/f+i. 
Then u(^r<i<i^i'^i) = by the recursive hypothesis. Besides, Prb~^u{Y2r<i<i ^i'^i) — 
because of Lemma Hence it(X]r<j<; HiVi) = which completes the recursive step. 

Summarizing, the following main theorem of the paper is proved. 

Theorem 4.9 Any (non-holonomic) module in {F[dx,dy]y of the differential type at least 1 
has a fractional- derivatives series non-zero solution. 

Remark 4.10 One could consider an ideal J C F[dxi, • • • ' ^xt] ^^^^^ of the differential type 
1 with a number of derivatives t > 2> and ask whether J has always a fractional- derivatives 
series solution? The answer to this question is negative already for t = 3 and an ideal J = 
{PI1P2) C F[dx^, dx2, dx^] (being generic of the differential type 1) generated by an operator 
pi of the first order and p2 of the second order. 

5 Duality between non-holonomic ideals and fractional- 
derivatives series solutions 

There is a well-known duality [13] between (left) differential ideals and their spaces of solutions 
(being an analogue of the duality between radical ideals and varieties in algebraic geometry) . 
To establish a similar duality for non-holonomic ideals in F[dx,dy] (so, of the differential 
type 1) we need to make use of the equivalence relation on ideals introduced in [10]. We 
say that non-holonomic ideals 7^ J, Jq C F[dx,dy] are equivalent if the leading coefficients 
of degree 1 (see Section Hilbert-Kolchin polynomials of three ideals J, Jq, J Ci Jq coincide 
(denote these leading coefficients by e), then moreover, GCD{symb{J)) = GC D{symb{Jo)) = 
GCD{.symb{JriJo)) and the degree of the latter polynomial equals e (see Lemma ITT]) . Li this 
case ideal (J, Jo) is also non-holonomic and GCD{symb{{J, Jq))) = GCD{symb{J)) as well 
[TOl O [21] and moreover, clearly four ideals J, Jq, J fl Jq, {J-, Jo) are equivalent. Equivalence 
classes of ideals play a similar role to classes (in algebraic geometry) of plane curves with the 
same sets of 1-dimensional components. Li this Section we prove that the sets of fractional- 
derivatives series solutions of equivalent non-holonomic ideals coincide and that there is a 
duality between the equivalence classes of ideals and their respective sets (which basically 
means that to distinct classes correspond distinct sets). 

In this Section we keep the notations from Section [H The next lemma states that the 
multiplication by 7^ p G F[dx,dy\ on D-module F = is an epimorphism (the conditions 
on its injectivity follow from Theorem 12.51 cf. also Remark 14. 4p . 
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Lemma 5.1 For any 7^ p G F[dx, dy] we have pV = V . 

Proof. Let Hi = hG^'o) + /^G^^n-*/^) g y (see We search for H = 

liQi^) -\- Y2i>i /ijG^*"*/''^ € V such that pH = Hi. Treating h, s as indeterminates we get from 
Lemma 11.41 that 

p(/iG(^)) =r/G(^+'^/'?) + ^.g(^+''/'?-'^/9) (7) 

l<j<K 

for certain integer k > and Unear ordinary differential operators t] ^ 0, rji in h (with 
coefficients being differential polynomials in fi,...,fk which we recall are assumed to be 
fixed). Indeed, such k with non-zero r/ exists since in the expansion of p{hG^'^^) the coefficient 
at G^*) in ([7j) equals p{h), so is a non-zero linear ordinary differential operator. Therefore, 
we put s = sq — k/q, and there exists /i G F for which r]{h) = h (because F is differentially 
closed). At the next step comparing the coefficients of pH and Hi at G^^~^'^^'^~^^'^^ one can 
find hi from an equation of the form r]{hi) = h for certain h F, and so on one can find 
/i2 , /i3 , . . . consecutively. ■ 

First consider two equivalent non-holonomic ideals J, Jo C F[dx,dy]. We claim that four 
sets of all fractional-derivatives series solutions of J, Jq, J n Jq, (J, Jo), respectively, coincide. 
As in Section [3] one can suppose w.l.o.g. that dyfi does not divide g = GGD{symb{J)) 
and consider left ideals J, Jq C being principal. Let J = (p), Jn Jq = (po) for suitable 
generators p,po G F[dx,dy]. Then bpo = p2P for appropriate p2 G F[dx,dy], 7^ 5 G F[dy]. 
Lemma 14.51 implies that symb{po) coincides with symb{p) up to a power of dyfi, whence 
symb{p2) is a power of dyfi. Therefore, Proposition 14.61 and Remark 14.41 entail the required 
claim on coincidence of the sets of fractional-derivatives series solutions of ideals J, J (1 Jq 
(and in a similar way also of Jq, (J, Jq))- 

Now let non-holonomic ideals J £ Ji C F[dx,dy] be non-equivalent. Our purpose is 
to find a fractional-derivatives series solution of J being not a solution of Ji. Denote gi = 
GG D{symb{Ji)) , clearly gi\g. Again we suppose w.l.o.g. that dyfi does not divide g. Let J = 
(p), Ji = (pi) C R. Then bip = P3Pi for suitable ps G F[dx,dy], / 61 G F[dy]. Since J, Ji 
are not equivalent we have deg{gi) < deg{g) (see Lemma HTJ and because of that Lemma 113] 
implies that symb{p3) has a divisor of the form dxfi + a^dyfi for a certain 03 G F. Due 
to Theorem 12.51 there exists a fractional-derivatives series solution Hi = ^^y^^ hiG^~'^/'^^ (see 
([2])) of equation p^H = for appropriate G = Gs2,...,Sk{fi: f2, ■ ■ ■ , fk) where dxfi + a^dyfi = 0. 
Now we apply Lemma 15.11 to pi and obtain a fractional-derivatives series H G Vq such that 
piH = Hi. Therefore, bipH = and hence pH = in view of Lemma |4.3[ Thus, H is a 
desired solution of J being not a solution of Ji. 

Finally, consider non-equivalent non-holonomic ideals J, Ji C F[dx,dy] and assume that 
their respective sets of fractional-derivatives series solutions coincide. Then ideal (J, Ji) 
has also the same set of fractional-derivatives series solutions, in particular (J, Ji) is non- 
holonomic by virtue of Remark 14.21 and of Theorem 12.51 Therefore, due to the proved above 
three ideals J, (J, Ji), Ji are equivalent which contradicts to the assumption. 

We summarize the proved duality in the following 

Proposition 5.2 Non-holonomic left ideals J, Jq C F[dx,dy] are equivalent if and only if 
they have the same sets of all fractional- derivatives series solutions. 
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For a non-holonomic ideal / C Fldx^dy] denote by [/] the equivalence class of non- 
holonomic ideals which contains / and by V{I) the set of all fractional-derivatives series 
solutions of /. In [10] we define the following partial ordering on the classes: [J] is subordi- 
nated to [/] if there exist ideals Ji G [J] , I\ € [/] such that Ji C /i . 

Corollary 5.3 [J] is subordinated to [I] if and only ifV{I) C V{J). 

Proof. Let V{I) C V{J), then V{{I,J)) = V{I). Proposition O entails that (/, J) is 
equivalent to /, hence [J] is subordinated to [/]. 
The inverse implication is evident. ■ 

Now we connect the subordination relation with localizations of ideals in the ring R = 
F[dx,dy\{F[dy\)-^ (see SectionS]). 

Proposition 5.4 // V{I) C V{J) then J <Z I d R (provided that dyfi does not divide 
GCD{symb{I))). 

Proof. Let J = {p), I = {q) for suitable p,q & F[dx,dy] (cf. above). Then p = poq for 
appropriate po G R- Whence V{I) C V{J) relying on Proposition 14.61 ■ 

For a pair of left ideals J C Ji C F[dx-i^, . . . , dx^] we have introduced in [10] a concept 
of relative syzygies. Namely, let Ji = {pi, . . . ,pt), then we define the left module of relative 
syzygies 

Syz{J,Ji) = {{qi, ...,qt): ^ qiPi e J; qi e F[dx^, . . .,dxj, l<i<t}. 

l<i<t 

Making use of [20j one can verify [10] that module Syz{J,Ji) is independent of a choice of 
generators pi, . . . ,pt. Let us denote by U{J) C F the space of solutions of J which can be 
treated as a C-vector space. It was proved in [10] that the quotient U{J)/U{Ji) is isomorphic 
to U{Syz{J,Ji)) C FK 

Here we establish a similar result for non-holonomic ideals J C Ji C F[dx,dy] and their 
spaces of fractional-derivatives series solutions VciJ) C Vg of the form ([2]) for any G fixed 
for the time being (see Section Hj), again we treat Vg'(^) as a C-vector space. As in Section H] 
one can assume w.l.o.g. that dyfi ^ 0. 

Mapping i\) : v ^ {pi, . . . ,pt)'^v assures a monomorphism Vg{J)/Vc{Ji) ^ 
VG{Syz{J, Ji)). To show that it is an epimorphism take an arbitrary vector (wi, . . . , wt) G 
VG{Syz{J, Ji)) C Vq. The following property holds: for any qi, . . . ,qt € F[dx, dy] such that 

Yli<i<t1iPi = (moreover, one can suppose that Yji<i<t1iPi ^ J) ^^'^^ Y^i<i<t H'^i = ^■ 
Clearly, this property holds also for any qi . . . ,qt £ R (see Section H]). Consider principal ideal 
Ji = {pi,... ,pt) = {p) C R for suitable p G F[dx,dy]. Then p = b^^ Z]i<i<t for appropri- 
ate / 6 G F[dy], rji,. . . ,r]t G F[dx,dy]. Denote w = J2i<i<t^i'^i ^ ^G- Due to Lemma [5!T] 
there exists v £ Vq such that (X]i<i<t Wi)^ = For each 1 < io < t one can find Ajg G R 
for which pi^ = Xi^p, the mentioned above property implies that Wi^ = \o^~^('l2i<i<t'ni'^i)i 
hence vJiq = PiqV, i. e. ij{v) = {wi, . . . ,wt). Finally, we check that v G Vg{J)- Indeed, an 
arbitrary q £ J can be represented as q = Ylii<i<t1iPi fo'^ certain qi, . . . ,qt G F[dx,dy], then 
{qi,...,qt) G Syz{J, Ji), whence qv = '}Zi<i<t1i'^i ~ ^- Thus, in the introduced notations 
we have proved the following 
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Proposition 5.5 For any non-holonomic ideals J C Ji C F[dx,dy] and G there is an iso- 
morphism of C-vector spaces Vg[J) /Vg{Ji) and VG{Syz{J, Ji)). 

One can deduce Proposition 15.21 from the latter Proposition invoking Theorem 14.91 
It would be interesting to clarify, whether for non-holonomic ideals /, J C F[dx,dy] the 
equality 

GCD{symb{I)) ■ GCD{symb{J)) = GCD{symb{{I, J))) • GCD{symb{I D J)) 

holds? Observe that the degrees of the polynomials in both sides of the latter equality coincide 
in view of [3l [2l] taking into account Lemma 14.11 A more subtle question is whether for any 
G the equality Vg{I n J) = Vg{I) + Vg{J) is true? 

6 Completeness of fractional- derivatives solutions for separa- 
ble linear partial differential operators 

Let T = Tn + • • • + Tq € F[dx,dy] be a separable LPDO, i.e. its symbol symhiT) = Tn = 
Y\i<i<nidxf — CLidyf) is the product of n pairwise distinct homogeneous linear forms in 
dxfjdyf. One can always bring symb[T) to this form monic with respect to dxf making, 
if necessary, a C-linear transformation of dx-, dy in case when symb{T) has a divisor dyf . 

For each 1 < i < n the equation T = has a fractional-derivatives series solution of the 
form (due to Theorem 12.50 

^*G(0)(/,)+/imG('') (/«) + •■■ (8) 
where dxfi — didyfi = and h = h^^i satisfies the first-order LPDE 

{dxh - aidyh) +T;ri{Wh = (9) 

We observe that hj,i;j = 1, 2, . . . satisfy similar to Q equations with the highest (first-order) 
form J"-t\ /.) (dxhj^i — aidyhj^i), being not necessary homogeneous. 

From now on throughout this section we assume that F is the field of meromorphic 
functions in a certain domain M C C^, thus the coefficients of T belong to F. For a suitable 
point (xo, yo) £ M the series ([8]) can be rewritten as a formal power series in x— xq, y — yo- Our 
goal is to find a point (xq, yo) and look for solutions of T = as power series in x — xq, y — yo- 

We choose a point (xq, yo) G M such that all the coefficients of T at this point are defined 
and in addition, the values ai(xo,yo) are pairwise distinct for 1 < i < n. The latter is 
equivalent to that the discriminant of symb{T) does not vanish at this point. Therefore, all 
the points of M out of an appropriate analytic subvariety of M of the dimension 1 satisfy 
these requirements. 

One takes a solution /j (being a power series in x— xq, y— yo) of the equation dxfi—aidyfi = 
with a vanishing free coefficient (which we denote by fi{xQ,yo) = 0) and with a non- 
vanishing vector of coefficients at the first powers of x — xo,y — yo (which we denote by 
{dxfi,dyfi){xo,yo), thereby dy/i(xo,yo) 7^ 0). We observe that this LPDE has always a 
solution with arbitrary chosen free coefficient and non- vanishing vector of the coefficients at 
the first powers of x — xo,y — yo since the vector of the coefficients (1, — Oj) at its highest 
(first) derivatives does not vanish at the point (xq, yo)- Hence the free coefficient of the power 
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series dxfi — o-jdyfi does not vanish when j ^ i due to the requirement on the discriminant. 
Therefore, by the same token one can find a solution h of the equation Q with a non-zero 
free coefficient which we denote by h{xo,yo) 7^ 0. 

We take an arbitrary solution of T = 0, being a power series in x — xq, y — yo and intend 
to represent it as a sum of n solutions of the form ([8]) (for 1 < i < n) in which IS 
replaced by its specialization (see Remark 1 1. 3 p 

j>o ■'' 

with indeterminate coefficients Cj^j G C. Then 

j>0 ^-^ ^ > 

Suppose that by recursion on k the coefficients Cj^i for j < k — 1,1 < i < n are already 
produced. Our purpose is to produce c^^j, 1 < i < n. Clearly, any solution of T = being 
a power series of the form Yip q>QW,q{^ ~ ^oYiv — VoY is determined by the coefficients hp^q 
with < p < n — 1. 

For each < p < n — 1 the contribution of the term at Ck,i (see ([8])) into hp k-p equals to 

h{xo,yo){dldy'P-^){xo,yo) = h{xo,yo){a'^{dyfi)''){xo,yo) 

taking into account that fi{xo,yo) = 0,dyfi{xo,yo) / 0, /i(xo,yo) / 0. 

Therefore, we obtain a linear (algebraic) system (in general, not necessary homogeneous) 
on Cfc^j;0 < i < n — 1 with the matrix being of the van-der-Monde type (a^(xo,yo))- This 
allows one to find uniquely Ck^i', < i < n — 1 and thereby, carry out the recursive step. 

Theorem 6.1 For a separable LPDOT of the order n with the coefficients being meromorphic 
in a certain complex domain M the sum of n spaces of specialisations of fractional- derivatives 
series solutions ofT = of the form (for fixed fi and hj^i) coincides with the space of all 
the solutions 0/ T = as formal power series in x — Xo,y — yo for any point (xq, yo) from M 
out of a suitable analytic subvariety of the dimension 1. 

It would be interesting to extend this theorem to a non-separable LPDO. Let us also 
mention that in [9j an algorithm for factoring a separable LPDO was produced. 

7 Applications to studying first-order factors of a linear par- 
tial differential operator 

7.1 Finding first-order factors of a linear partial differential operator 

Let T = Tn + ■ ■ ■ + Tq he Su LPDO of an order n in 2 independent variables, where Tj = 
Yi'^i,j-idxdy~^ is a sum of the derivatives of the order j. We assume that the coefficients 
cii^j are taken from the field Q(x, y) in order to design algorithms, while / is taken from a 
universal field F (cf . Section [6]) . 
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As we are looking for the first-order factors of T of the form L = dx + ady + b € F[dx, dy\ 
we need to study the solutions of L = (w.l.o.g. one can assume that the coefficient at dx of 
L does not vanish, otherwise one can change the roles of x and y). Take any solution / of the 
symbol {dx + ady)f = of L such that dyf / and consider G = G^'^\f) (cf. Section E]). For 
any h G F being a "particular" solution of L = 0, we have that hG is a fractional-derivatives 
series solution of L = 0. 

Lemma 7.1 An operator T has a right first-order factor L if and only if the equation T = 
has a solution of the form hG. 

Proof. If T has a right factor L then T has a solution hG. 

Conversely, assume that T = has a solution hG. Dividing T with remainder by L one 
can represent T = SL + J2o<i<n^i^y ^ suitable operator S. Consider the largest k such 
that bk 7^ 0. Then in the expansion of {Eo<i<kbidi)hG in {G^'^ the coefficient at G^'^) 
equals bkh{dyf)^ ^ 0. The obtained contradiction shows that T = SL. ■ 

Thus, we are looking for a solution of T = of the form hG. Expanding T{hG) = AqG^^^ + 
■ ■ ■ + AnG^'^\ we get first that An = symb{T). Therefore, we fix for the time being a linear 
divisor of the form dxf + adyf of symb{T) and assume that this divisor vanishes. Thereby, 
the calculations below (arthmetic manipulations and polynomial factoring) will be carried 
out over the field Q{x,y)[a]. This can be fulfilled representing Q{x,y)[a] ~ Q{x,y)[z]/{g) 
where g G Q{x,y)[z] is the minimal polynomial of a (see [5]). So, we obtain n equations 
Aq = ■ ■ ■ = An-i = treated as LPDO in h with the coefficients being non-linear differential 
polynomials in /. We denote the ring of all these polynomials by P = Q{x,y)[a]{dxf,dyf}. 
Applying to Aq = ■ ■ ■ = A^-i = the procedure of constructing a Janet base [22j one gets the 
conditions of solvability in h of Aq = ■ ■ ■ = An-i = expressed as a disjunction of systems of 
the form 

Pi = ••• =M = 0,po ^0 (10) 

where pi € P. Using the relation dxf + adyf = one can reduce each pi to an (ordinary) 
differential polynomial pi in dyf. Denote the ring of ordinary differential polynomials by 
R = Q{x,y)[a]{dyf}. 

Applying to the formula pi = ■ ■ ■ = pi = 0,po 7^ the subroutine of the elimination 
procedure in the theory of ordinary differentially closed fields from [23] (see also [6] where 
its improvement with a better complexity bound was designed) one obtains an equivalent 
disjunction of systems of the form 

r = 0,ro/0 (11) 

for suitable differential polynomials r, rg G R. Briefly, this subroutine consists in alterna- 
tive executing 2 types of steps while there are more than one equality of (ordinary) dif- 
ferential polynomials. The first type of steps is executed when all the highest derivatives 
occurring in these polynomials are equal, in this case the algorithm calculates their CCD 
viewing them as (algebraic) polynomials in this highest derivative (and branching depending 
on vanishing the leading coefficients). Else, if not all the highest derivatives are equal, as 
the second type of steps one can diminish the highest derivative. Moreover, if r contains 
the dyf as its highest derivative then r considered as an (algebraic) polynomial in the ring 
K = Q(x, y)[a] [/, dyf, . . . , dyf] is irreducible. In addition, tq is less than r with respect to the 
term ordering, i. e. if contains dy° f as its highest derivative then either ko < k or ko = k 
and the degree of rg with respect to dyf is less than the similar degree of r. 
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Replace dxf by —adyf in d^r. This yields a differential polynomial f G i? of the order 
at most k + 1 (its role is similar to an 5'-pair in Janet type algorithm [22] )• If r does not 
belong to the differential ideal (r) C R, we again apply to the system r = f = 0, rg 7^ the 
used above subroutine from the elimination procedure and get an equivalent disjunction of 
systems of the form (lllh with less term ordering than of r and continue as above. 

Now assume that r belongs to (r). Then we claim that any solution of (jlip provides a 
solution of (fTOl) . Indeed, otherwise, the ideal (r, dxf+adyf) C P would contain an appropriate 
power Tq [13], p. 146-148. This yields a relation of the form 

ro = Y^ Aijd^xdi^ + Bi,jdidi{dxf + adyf) 

for suitable Aij,Bij G P. Replacing in this relation dxf for —adyf and taking into account 
that f belongs to (r), we deduce that 

r^ = 5]i,4r (12) 
j 

for certain Aj G R. From the equation dyr = we express 



d{d^J) 



for an appropriate Bk+i G K. After that express successively 



didifY y ' d{dify-- 

Substitute these expressions in this results in the equality 

for some t and A E K. But r is irreducible in K and ro is less than r with respect to the 
term ordering. The obtained contradiction proves the claim and the following theorem. 

Theorem 7.2 There is an algorithm which tests whether an operator T G Q{x, y)[dx, dy] has 
a first-order factor with the coefficients in a universal field F. The algorithm invokes two 
subroutines: the elimination of an unknown function in a system of LPDO's (in other words, 
a parametric Janet base), and a subroutine from the elimination procedure in the theory of 
ordinary differentially closed fields. 

Remark 7.3 // one uses a direct method of finding the coefficients of a first-order operator 
L and of an (n — l)-th order Q such that T = QL, then one has to apply an elimination in 
the theory of partial differentially closed fields whose complexity is unclear how to estimate in 
a reasonable way (cf. f^]. 

Remark 7.4 One can also search for left first-order factors of an LPDO (by means of con- 
sidering an adjoint operator). 

Corollary 7.5 There is an algorithm to factor LPDO's of the orders at most 3. 
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7.2 Intersection of principal first-order ideals 

In this subsection by F we denote a differential field with derivatives dx,dy. 

First consider the ideals li = {dx + ady + bi) with the same highest (first-order) forms 
where a,bi € F,l < i < n. 

Proposition 7.6 The ideal Ji n • • • H /„ is principal 

Proof. Denote E = d^ + ady. The ring F[E] is left-euclidean, therefore, the intersection 
ii n . . . n In = (Q) C F[E] is principal where we denote Ij = {dx + ady + hi) C F[E] and 
Q = QsE^ + • • • + go for certain qq, . . . ,qs & F,qs ^ and s < n. 

Our aim is to prove by induction on n that /i n ■ • ■ H I„ = {Q) C F[dx, dy]. Assume that 
it is already proved and consider the intersection /i n • • • fl /„ H /„+i. There can occur two 
cases. Either /i n . . . fl H In+i = /i fl . . . fl in this case Q = N{dx + ady + hn+i) for a 
suitable N G F[E], therefore, D (Q) and /i n • • • n /„ n = (Q). 

Or else /i n . . . n /„ ^ /"i n . . . n /„ n = (M) for an appropriate M = ms+iE''+^ + 
■ ■ ■ + mo € -^[-S] with mo, • • • , m^+i G F. Clearly, M G /i fi • • • fl n /„+!. It is necessary 
to show that for any F G Ii n • • • fl /„ fl In+i we have F G (M). Since the highest derivative 
with respect to dx which occurs in M is d^.^^, one can divide V hy M with remainder and 
get V = WM + U where W,U e F[dx,dy] for a certain [/ G /i n • • • n /„ n such that 
So = ordd^{U) < s. If [/ = we are done, so suppose that U ^ 0. We have 

U = ZQ = T{dx + ady + bn+i) (13) 

for suitable Z,T & F[dx,dy], hence sq = s and ordd^{Z) = 0,ord(i^{T) = s — 1. One can 
expand T = ts-iE^^^ + • • • + fo for appropriate to, ... , tg-i G F[dy]. Thus, the equation (fT3]) 
one rewrite with respect to the powers of E: 

Z{qsE' + . . . + go) = {ts^iE'-^ + --- + to){E + bn+l) 

which is equivalent to a system of the following s + 1 equalities: 

~ + ^j^3,j + + ■ ■ ■ + is-i^i,s-i (14) 

for suitable 6^^, . . . , G F; 1 < j < s and 

Zgo = tobn+l + tl6o,l + ^2^0,2 H h ts-l&0,s-l (15) 

Viewing the right-hand sides of the equations ()14p . ()15p as a linear system in io; • • • > ^s-i we 
get that there is a unique linear combination (from the right) of s expressions in the right- 
hand sides of (fT^ which equals the coefficients /i, . . . , of this combination belong to 
F. Therefore, the solvability of (fHI) . (fTSll in Z 0,to, . . . ,ts-i entails the equality 

qifi H h g^/s = qo (16) 

Thus (jl3p implies (I16p . Hence as a solution of the system (I14p . (jlSp one can take Z = 1 
and consecutively express t^-i ^ F from the equation with j = s, after that express 
ts-2 £ -P" from the equation (I14p with j = s — 1 and so on, finally express to & F from 
with j = 1. The last equation (fT5]) of the system is fulfilled due to (fT6|) . As a result 
we obtain (cf. ([13])) Q = (^s-i-E'"'"^ + ••• + to){E + bn+i) with G F, in other words 

iin...n/„ = (g) c/„+i cF[F]. 
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This leads to contradiction with the assumption Ii n . . . fl /„ ^ Ii D . . . f] In In+i, which 
shows that the supposition U ^ was wrong, thus Ji fl . . . n/„ n/„+i = (M). The proposition 
is proved. ■ 

Corollary 7.7 The ideal /i H . . . fl /„ is generated by an element from F[E]. 

Now let the ideals Ij = {dx + aidy + 6j) C F[dx,dy] be given, where ai,bi & F,l < i < k. 
Our goal is to study their intersection I = IiD . . . Ik- Combining together all the classes of 
the ideals with the same Oj and making use of Corollarv 17.71 we replace the intersection from 
one class by (Zi) for a certain Z,, E F[Ei] where Ei = d^ + a^dy. Then / = /i fi . . . n = 
{Zi) n • • • n {Zi) for some I. Denote Si = ord{Zi); 1 < i < I and s = si + • • • + s/. 

Lemma 7.8 For any Q ^ I we have ordd^{Q) > s. 

Proof. Observe that symb{Q) is divided by ni<j</(^^/ + o-idyfY^ treated as a homoge- 
neous polynomial in dxf,dyf. ■ 

Theorem 7.9 a) The ideal I is principal if and only if I contains Q with the order ord{Q) < 

s; 

b) in this case ord{Q) = s and I = {Q). 

Proof. Obviously, the typical differential dimension dim{{Zi)) = Si;l < i < I [13j and 
dim{I) < s due to [3l|2l]. Hence if / = (L) is principal then ord{L) = dim{I) < s. 

Conversely, let Q € / and ord{Q) < s, by virtue of Lemma 17.81 we have ord{Q) = s and 
the derivative occurs in Q. Our purpose is to show that / = (Q). Indeed, take any V & I 
and divide V hy Q with remainder, we get V = WQ + U where ordd^{U) < s, therefore, 
U = due to Lemma 17.81 Thus, I = (Q). ■ 

Corollary 7.10 let the differential field F = Q{x,y). There is a polynomial-time algorithm 
which tests whether I is principal. 

Proof. First the algorithm produces Zi;l < i < Ihy finding a non-zero solution of a linear 
(algebraic) homogeneous system on the coefficients from F of Ti, . . . ,Tn € -Fi-Ej] such that 
Ti{dx + Oidy + 6i) = ■ • • = Tn{dx + aidy + bn) with the minimal possible order ord(Ti) = ■ • • = 
ord(Tn) (trying consecutively the orders 1,2...). Denote Zi = Ti{dx + Oidy + 6i), Si = ord{Zi), 
then Zi is a generator of the ideal {dx + Oidy + 6i) n • • • fl {dx + aidy + see Corollary 17.71 

Thereupon the algorithm looks for 14, . . . , S F[dx, dy] with ord{Vi) < s — Si;l < i < I 
such that V\Z\ = ■ ■ ■ = ViZi. The latter we treat as a linear (algebraic) homogeneous system 
in the coefficients from F of Fi, . . . , VJ. Theorem 17.91 entails that this system has a non-zero 
solution if and only if / is principal. ■ 

Remark 7.11 Observe that the usual method of finding the intersection of ideals invoking 
Groebner bases, runs in double- exponential time. 
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7.3 Constructing intersection of all first-order factors 



In this subsection F denotes a universal field [13] with two derivatives dx,dy. 

The purpose of this subsection is to construct the intersection U C F[dx,dy] of all the 
principal ideals (L) for the first-order factors L G F[dx,dy] of T € Q{x,y)[dx,dy]. Evidently, 
U D (T). We mention that in [TO] a radical of a module of a differential type r was defined 
as the intersection of the maximal classes of r-equivalent modules, and a question was posed 
whether one can calculate the radical. Here U (which could be called a first-order radical) is 
defined as an ideal (rather than a class of equivalent ideals) and moreover, we calculate U. 

Observe that the construction from the Subsection 17.11 represents the family V of all the 
solutions of the form hG (and which correspond to first-order factors of T due to Lemma [7. ip 
as follows (we use the notations from Subsections 17. H 17.20 . We assume that a is fixed, while 
/ just satisfies the equality dxf + adyf = 0. The family y is a union of subfamilies of the 
form Vq where Vq is given by means of a Janet base 



for h where Vi^^i^^i G R together with a system ([TT]) for /. 

For each element hG G V consider the first-order LPDO Lhc = dx + ady + bhc such 
that LfiQ{hG) = (see Lemma 17. 1|) . We claim that one can extend Proposition 17.61 from 
a finite to an infinite number of principal ideals and conclude that the ideal rihGeviLhc) is 
principal and moreover, is generated by a suitable element Q = J2o<i<s li^^ ^ ^[^] (see 
Corollary I7.7p . Indeed, one add consecutively the ideals Ii = {Lh^Gi)-,!^ = {Lh2G2)-> ■ ■ ■ for 
hjGj G V , while the intersection /i n • • • fl fl /j C n • • • n decreases (cf. the proof 
of Proposition 17. 6p . Then /i n ■ ■ ■ n /j = {Qj = Yl,o<i<j Qi,j^^) foi' appropriate qij G F (cf. 
the proof of Proposition 17. 6p . Hence (T) C h (1 ■ ■ ■ H Ij = (Qj) due to Corollary 17. 7[ Thus, 
j < n and nhGeviLhc) = h H ■ ■ ■ (1 Ij which proves the claim. 

To produce Q = Qj = J2o<i<j the algorithm successively tries j = 0, 1, . . ., treating 
Qi as indeterminates. The aim is to find Q such that Q{hG) = for any hG G Vq (for each 
subfamily Vq of V). The algorithm expands Q{hG) = AqG^^'^ + ■■■ + AjG^^'^ (cf. Subsec- 
tion [77T])- One can view each Ai as an LPDO in h with the coefficients being linear forms 
in qo,. . . ,qj over R. The algorithm divides every Ai,0 < i < j with the remainder by the 
Janet base (fT7|) . as a result we obtain LPDO Ai = Y^-_^ ai^i^^i^d^^d^^ . Thus, Q vanishes at 
any hG G Vq if and only if ai^i^^i^ = for all < i < j; ii, ^2 under condition (jlip . 

Denote by S the conjunction of the systems ai^i^^i^ = for all < i < j;ii,i2 and for all 
subfamilies of the form Vq of ^ . One can treat <S as a homogeneous linear over qQ, . . . ,qj system 
with parameters being derivatives /, dyf, . . . , dyf for a certain /. Solving this parametric linear 
system (see e.g. [7]) the algorithm finds the (algebraic) conditions on f,dyf, . . . ,dyf under 
which the system is solvable and in addition, finds the expressions for solutions (being rational 
functions in the parameters). After that the algorithm tests whether these conditions are 
compatible with (jlip . applying the subroutine from the elimination procedure which yields 
formula ([TT]) in Subsection 17.11 If yes then the algorithm produces a solution qo, . . . ,qj G R 
of the parametric linear system. Else, the algorithm proceeds from the current value j to the 
next value j + 1- 

Thus, the algorithm for each a such that dxf + adyf is a (linear) divisor of symb{T), pro- 
duces applying the described above construction a generator Qa G F[dx, dy] of the (principal) 
ideal being the intersection of all the principal ideals generated by the divisors of the form 




(17) 
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dx + ady + 6 of T for varying h. Finally, the algorithm finds the intersection U = r\a{Qa) 
over all the divisors dxf + adyf of symb(T) making use of Janet bases (cf. ilQj). Thus, the 
following theorem is proved. 

Theorem 7.12 For any LPDO T € Q{x , y)[dx , dy] one can construct the intersection of all 
the principal ideals generated by the first- order factors of T. 

8 Fractional-derivatives series solutions of a second-order op- 
erator and factoring 

In this section we study a particular case of a second-order LPDO T = Tq + Ti + T2 and 
describe its possible fractional-derivatives series solutions being outputs of the algorithm 
from Section El First, if the symbol symb{T) is separable then for each of its two different 
linear divisors dxfi + adyfi the algorithm provides a fractional-derivatives series solution of 
T of the form (cf. ([HD) 

0<j<oo 

where G = G{fi) and dxfi + adyfi = 0, /i 7^ const. 

From now on let us assume that symb{T) is non-separable and write T = d"^ + 2adxdy + 
a^dy + bo^idx + bifidy + 60,0 • The first step of the algorithm from Section [2] yields /i such that 
dxfi + o-dyfi = 0. Introduce the discriminant of T as follows: 

{-T + 6o,o)/i = {dxa + 2adya + a&o,i - bi,o){dyfi) = Disc ■ (dyfi). 

If Disc 7^ we take any /2 which satisfies the following (non-homogeneous) first-order 
LPDE: 

dxf2 + adyf2 = ^JDisc■ {dyfi) 

and G = Gi/2(/i,/2) (see Definition II. ip . then the algorithm constructs a fractional- 
derivatives series 

0<i<oo 

being a solution of T = 0. Each of two values of the sign of the square root provides a generic 
solution of the multiplicity 1 (see Section [3]). It corresponds to the leading edge with the 
endpoints (0, 2), (1, 0) having the slope 1/2 at the second step of the algorithm. 

When Disc = the algorithm yields a (fractional-derivatives series) solution hG{fi) of 
T = for an arbitrary particular h such that T{h) = 0. It corresponds to the leading edge 
with the endpoints (0,2), (0,0) having the slope at the second step of the algorithm and 
provides a generic solution of the multiplicity 2. Relying on Lemma 17.11 one obtains the 
following corollary (cf. [9J). 

Corollary 8.1 A second-order LPDO with a non-separable symbol is irreducible if and only 
if Disc / 0. 
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